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METHODS AND CIRCUITS FOR IMPLEMENTING 
PROGRAMMABLE GAMMA CORRECTION 

FIELD OF INVENTION 

[0001] The present invention relates in general to video processing and in 
particular to methods and circuits for implementing programmable gamma correction. 

BACKGROUND OF INVENTION 

[0002] Gamma correction is a well-known technique for controlling the luminance 
(brightness) of the display pixels produced by cathode ray tube (CRT) displays and the 
like. Generally, the actual luminance of each pixel on the display screen has a non- 
linear relationship with the con-esponding pixel data, such as the digital pixel data output 
from the DVD decoder. For most CRT displays, this non-linear relationship is the power 
function x\ in which x is the brightness of the given pixel defined by the pixel data, y is 
the gamma value for the given display system, and x^ is the actual brightness of the 
pixel produced by the display system. Normally, the values of x correspond to analog 
display control voltages between zero (0) and one (1) volt, and therefore for a given 
value of X, the brightness of the con-esponding pixel on the screen is typically too low. 
Gamma correction essentially corrects for the non-linearity in the conversion of the input 
pixel data into display pixels by applying the inverse function to the input data (i.e. x^'^ ) 
in either software, hardware, or both hardware and software. As a result, the overall 
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brightness of the display pixel, as well as the ratio between color components of that 
display pixel, is con-ected relative to the input data. 

[0003] The value of y depends on the particular display system, as well as the 
lighting environment surrounding that display system during operation. Typically, the 
value of y is between 1.4 and 3.5. The typical value of y in the case of conventional 
display is 2.5. for an Apple Macintosh display, y is taken to be 1 .8, while the NTSC 
standard is based on a y value of 2.2. Other systems, such as PAL and SECAM, 
assume a y value of 2.8. High quality video products (encoders/decoders) often also 
require that the error between the output of the gamma correction circuit and the ideal 
output, be less than two (2) least significant bits (LSBs). 
[0004] One widely-used gamma correction technique utilizes a look-up-table 
(LUT) in memory. In this case, gamma conrected values for all possible pixel data 
values from the corresponding data source are pre-computed and stored in the LUT. 
Each pixel value input from the data source indexes a gamma corrected output value in 
the LUT, which is then utilized to generate the corresponding pixel on the display 
system's display screen. The LUT approach, however, consumes a significant amount 
of memory space. Furthemnore, the LUT must be re-programmed each time a different 
gamma value is required for driving a different display system. 
[0005] In another LUT - based approach, difference values are calculated 
between the gamma corrected output values to the display system and the linear input 
values from the data source. These difference values are then stored in the LUT 
instead of the absolute gamma corrected output values discussed above. Since the 
dynamic range of the difference values is less than the dynamic range of the absolute 
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values, the size of the LUT table Is reduced. However, this second approach is still not 
memory efficient or easily programmable. 

[0006] In a further gamma correction technique, a curve representing the gamma 
correction operation is approximated as piecewise linear line segments. The location of 
the knees or edges of each piecewise segment is fixed and the output value at the 
knee-points are programmed in a register bank. This approach introduces significant 
enror, especially for lower luminescence values and also requires the re-programming of 
registers with pre-computed values to update the gamma correction curve. 
[0007] Similarly, in one additional technique, a difference curve is approximated 
between the gamma corrected output values and the linear input values as piece-wise 
linear segments. The corresponding segment information including slope, intercept, 
and segment end points (knee-points) are stored in a LUT or register bank. This 
technique has improved error characteristics but still requires a significant amount of 
programming, after pre-computing the segment information, to change the gamma 
value. 

[0008] In sum, new gamma correction circuits and methods are required, which 
support gamma correction in a wide range of systems and which provide ease of 
programmability, while at the same time minimize the introduction of error. 



SUMMARY OF INVENTION 

[0009] A method of gamma correction includes selecting lower and upper 
reference curves corresponding to selected gamma values. A gamma correction curve 
is generated from a corresponding gamma correction value and cross-correlated with 
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the upper and lower reference curves to generate a corresponding set of cross- 
correlation factors. The set of cross-correlation factors are stored and indexed to the 
corresponding gamma value. An input value is received for gamma correction with the 
corresponding gamma value. Data from the upper and lower reference curves indexed 
by the input value is then operated one with the cross-correlation factors to generate a 
gamma corrected output value. 

[0010] Embodiments of the principles of the present advantages afford 
substantial advantages over the prior art. The use of reference curves and cross- 
correlation factors allow a range of gamma values to be programmed with minimal 
memory space. Specifically, the only memory spaces required are those necessary for 
storing the reference curves and the cross-correlation factors; memory space is not 
required for storing a large set of complete gamma correction curves since a range of 
gamma correction values are implicit in the stored set of cross-correlation values. 
Further, the desired gamma value for a given application can be directly selected, 
simply by writing to a single register. Additionally, the minimum error between the 
actual output values generated in accordance with the inventive principles and the 
corresponding ideal output values is reduced in view of conventional gamma correction 
schemes. In particular, the correlation-based (i.e. non-linear-based) interpolation of the 
present invention provides a substantially lower minimum error versus an ideal gamma 
correction curve than does conventional linear interpolation. 
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BRIEF DESCRIPTION OF DRAWINGS 

[001 1] For a more complete understanding of the present invention, and tiie 
advantages tliereof, reference is now made to the following descriptions taken in 
conjunction with the accompanying drawings, in which: 

FIGURE 1 is a high level block diagram of a video system suitable for 
describing one application of the principles of the present invention; 

FIGURE 2 is a flow chart illustrating an exemplary procedure for 
generating one or more desired gamma correction curves according to the present 
inventive principles; and 

FIGURE 3 is a block diagram of a programmable gamma correction circuit 
embodying the principles of the present invention and suitable for such applications as 
the gamma connection block of FIGURE 1. 



WSM Docket No. 
2836- P277US 



ATTORNEY DOCKET NO. 6 PATENT 

1434-CA 



DETAILED DESCRIPTION OF THE INVENTION 

[0012] The principles of the present invention and their advantages are best 
understood by referring to the illustrated embodiment depicted in FIGURE 1- 3 of the 
drawings, in which like numbers designate like parts. 

[0013] FIGURE 1 is a high-level block diagram of a video system 100 such as a 
DVD player suitable for describing one application of the principles of the present 
invention. Generally, DVD player 100 includes a digital data source, such as DVD 
processor 101, a digital video encoder 102, and a display system 103, such as a 
cathode ray tube (CRT) display system. 

[0014] Encoder 102 inputs data and control signals from DVD processor 101, 
including digital pixel data, in a format such as the standard YCbCr and YUV formats, 
along with control signals including the vertical and horizontal synchronization signals. 
An input demultiplexer (demux) 104 splits the pixel luminance ("luma") (Y) components, 
the pixel chrominance ("chroma") (CbCr or UV) components, and the timing and control 
signals. 

[0015] The luma components are processed through a data path including 
gamma correction circuitry 105, embodying the principles of the present invention and 
discussed in further detail below, and luminance ("luma") processing block 106. 
Luminance processing block 106 performs such operations as interpolation, sharpness 
enhancement, brightness, and contrast adjustments. Similarly, the chroma components 
are processed through parallel chroma processing path 107, which performs operations 
such as interpolation, filtering, and modulation. Timing generator block 108 processes 
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horizontal and vertical synchronization timing to generate a composite synchronization 
signal. 

[0016] The outputs from luma processing path 106, chroma processing path 107, 
and timing generator block 108 are multiplexed by output multiplexer ("mux") 109 into 
the inputs of output digital-to-analog converters (DACs) 1 10. DACs 1 1 0 provide the 
analog signals required to drive display system 103. Depending on the particular 
configuration of display system 103, DACs 1 10 output analog data in a format such as 
composite video, S-video, or RGB video. 

[0017] Programming interface and registers block 1 1 1 supports the 
programmability of various operations of encoder 102 by an external controller. For 
purposes of the present discussion, programming interface and registers block 1 1 1 
allows an external controller to program gamma correction block 105 in accordance with 
the inventive principles discussed below in conjunction with FIGURES 2 and 3. 
[0018] FIGURE 2 is a flow chart illustrating a procedure 200, which generates 
one or more sets of cross-correlation matrices representing one or more desired 
gamma correction curves, according to the present inventive principles. The cross- 
correlation matrices generated in accordance are then stored for later utilization during 
gamma correction of an incoming luma stream by interpolation between a pair of stored 
reference curves. One particular circuit suitable for implementing procedure 200 is 
discussed further below in conjunction with FIGURE 3. 
[0019] For the current cross-correlation matrix being generated, the 
corresponding gamma value G is selected at Block 201 . In the illustrated embodiment, 
G takes on values between 1 .4 and 3.5 in increments of 0.01 . 
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[0020] In accordance with procedure 200, at least one pair of precalculated and 
stored upper and lower reference curves are selected in response to the current gamma 
value G. In the illustrated embodiment, if G is less than or equal to 2.3 at Block 202, 
then at Block 203, upper and lower reference curves U(n) and L(n) are selected 
corresponding to reference gamma values of 2.3 and 1 .6 indexed by input values 
LUMA_iN(n). In addition, a desired gamma curve is generated from the desired gamma 
value G and input values LUMA_iN(n). The desired gamma curve is generated for 
calculating the cross-correlation matrices discussed below, although the desired gamma 
curve is not stored for utilization in actual gamma correction operations since the 
desired gamma curve is implicit in the corresponding cross-correlation matrices. 
Specifically: 



U(n) = gamma(LUMAjN(n), 2.3) = [LUMA_iN(n)]^'^-^ ; 
L(n) = gamma(LUMA_iN(n), 1 .6) = [LUMAjN(n)]^'^-^ ; and 
l(n) = gamma(LUMAjN(n), G) = [LUMA_iN(n)]^'°. 



[0021] On the other hand, in the illustrated embodiment, if at Block 202, the 
cunrent gamma value G is greater than 2.3, then a second set of pre-calculated and 
stored curves is selected at Block 204. In this example, the upper and lower reference 
curves U(n) and L(n) respectively correspond to reference gamma values of 3.3 and 
2.3. A curve representing the application of the desired gamma correction value G to 
the input values lumajn is generated for cross-correlation matrix calculation purposes. 
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In particular: 

U(n) = gamma(LUMAjN(n), 3.3) = [LUMAjN(n)]^'^-^ ; 

L(n) = gamma(LUMAjN(n), 2.3) = [LUMAjN(n)]^^-^ ; and 

l(n) = gamma(LUMA_iN(n), G) = [LUMA_iN(n)]^'°. 

[0022] At decision Block 205, for input values of LUMA_iN(n) less than 32, a first 
set of cross-correlation matrix entries are calculated at Block 206. The cross-congelation 
entries for the Illustrated embodiment are defined as follows. 
[0023] The matrix entry ^ILJ is the numerical value of the expected (mean) of 
the dot product of the vectors l(n) and L(n) from the curves defined at Block 203 or 204 
for values of LUMA_iN(n) less than 32 (i.e. n = 0,1 31 ). In particular: 

'^ILJ = E (l(n).L(n)) = Mean [l(0)*L(0) + l(1)*L(1) + ... + l(k). L(k)], 

in which k is a dummy variable between 0 and 31. 
[0024] The matrix entry *IU_I is the numerical value of the expected (mean) of 
the dot product of the vectors l(n) and U(n) from the curves defined at Block 203 or 204 
for values of LUMAjN(n) less than 32 (i.e. n = 0,1 ,...,31 ). In particular: 

*^IUJ = E (l(n).U(n)) = Mean [l(0)*U(0) + l(1)*U(1) + ... + l(k). U(k)], 
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in which !< Is a dummy variable between 0 and 31 . 

[0025] The matrix entry ''^ULJ is the numerical value of the expected (mean) of 
the dot product of the vectors U(n) and L(n) from the curves defined at Block 203 or 204 
for values of LUMAjN(n) less than 32 (I.e. n = 0,1 31 ) or: 

%L_\ = E(U(n).L(n)) = Mean[U(OrL(0) + ... + U(k).L(k)]. 

[0026] The matrix entry \l_I is the numerical value of the expected (mean) of 
the square of the vector L(n) from the curves defined at Block 203 or 204 for values of 
LUMAjN(n) less than 32 (i.e. n = 0,1 31): 

*^LLJ = E(L(n).L(n)) = Mean[L(0)*L(0) + ... + L(k).L(k)]. 
[0027] The matrix entry ^''UUJ Is the numerical value of the expected (mean) of 
the square of the vector U(n) from the curves defined at Block 203 or 204 for values of 
LUMA_iN(n) less than 32 (i.e. n = 0,1, ...,31). Specifically: 

*UU_I = E(U(n).U(n)) = Mean[U(0)*U(0) + ... + U(k).U(k)]. 

[0028] For values of LUMAjN(n) between 32 and 875 (i.e. n = 32 875), the 

following cross-correlation matrix entries are calculated at Block 207. 

[0029] The matrix entry *''lL_u is the numerical value of the expected (mean) of 

the dot product of the vectors l(n) and L(n) for the curves defined at Block 203 or 204 for 
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values of LUMAjN(n) between 32 and 875 (i.e. n = 32, 33,. ...875). In particular: 

"^IL.u = E(l(n).L(n)) = Mean[l(32rL(32) + l(33rL(33) + ... + IG).LO)], 

in which j is a dummy variable between 32 and 875. 
[0030] The matrix entry *''lU_u is the numerical value of the expected (mean) of 
the dot product of the vectors l(n) and U(n) for the curves defined at Block 203 or 204 
for values of LUMAjN(n) between 32 and 875 (i.e. n = 32, 33 875). In particular: 

*^IU_u = E(l(n).U(n)) = Mean[l(32)*U(32) + l(33)*U(33) + ... + IG).UG)], 

in which j is a dummy variable between 32 and 875. 
[003 1] The matrix entry *''UL_u is the numerical value of the expected (mean) of 
the dot product of the vectors U(n) and L(n) from the curves defined at Block 203 or 204 
for values of LUMA_iN(n) between 32 and 875 (i.e. n = 32,33 875) or: 

*^UL_u = E(U(n).L(n)) = Mean[U(32)*L(32) + U(33)*L(33) + ... + UG).LO)] 
in which j is a dummy variable between 32 and 875. 

[0032] The matrix entry ^LL_u is the numerical value of the expected (mean) of 
the square of the vector L(n) from the curve defined at Block 203 or 204 for values of 
LUMAjN(n) less between 32 and 875 (i.e. n = 32,33 875): 
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LL_u = E(L(n).L(n)) = Mean[L(32)*L(32) + L(33)*L(33) + ... + LO).L(j)]. 
in which j is a dummy variable between 32 and 875. 



[0033] The matrix entry **'UU_u is the numerical value of the expected (mean) of 
the square of the vector U(n) from the curve defined at Blocl< 203 or 204 for values of 
LUMAjN(n) between 32 and 875 (i.e. n = 32,33 875). Specifically: 



<D 



UU_u = E(U(n).U(n)) = Mean[U(32)*U(32) + U(33)*U(33) + ... + UG).U(i)]. 
in which j is a dummy variable between 32 and 875. 



[0034] The correlation matrix for values of lumajn below 32 is calculated, by 
utilizing the matrix entries calculated at Block 206 and the following matrix equation: 



aj = [*ILJ. '^lUJ* 



0 



LL I 



UL I 



UL I 



4> 



UU I 



in which aJ is a one by two matrix [a_l(1) aJ (2)]. 
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[0035] The correlation matrix for values of lumajn in the range of 32 to 875 is 
calculated at Block 209. using the matrix entries calculated at Block 207 from the matrix 
equation: 



r 



LL u 



UL u 



0 



UL_u 
UU u 



in which a_u is a one by two matrix [a_u(1) a_u (2)]. 
[0036] For each value of G programmed, four cross-correlation coefficients 
con-esponding to matrices a J and a_u are stored in memory of Block 210, for example 
in a LUT, for latter utilization during gamma con-ection operations on received values of 
LUMAJN. One particular exemplary LUT embodiment is discussed below in conjunction 
with FIGURE 3. 

[0037] During display generation, luma values are received at Block 21 1 for 
conversion. For values of lumajn less than thirty-two (32) at decision Block 212, a is 
taken as aj at Block 215, otherwise a is taken as a_u at Block 214. At Block 216 the 
gamma corrected value of a current input value lumajn by the desired gamma value G 
is calculated as: 



out = a* [L U]' 

in which [L U]' is the transpose of a matrix representing the lower (L) and 
upper (U) reference curves selected at Block 202. Alternatively expressed, each 
gamma corrected value for the desired gamma con-ection value G is calculated as: 
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out(n) = a(1)*L(n) + a(2)*U(n) 



[0038] in which ouT(n) is the gamma corrected value being calculated for the 
current input value LUMA_iN(n), L(n) is the value along the lower reference curve 
corresponding to the current lumajn(n), U(n) is the value along the upper reference 
curve conresponding to the current LUMA_iN(n), a(1) is the cross-conrelation matrix entry 
at the first column of the one by two (1 x 2) matrix calculated at Block 208 or Block 209, 
and stored at Block 210, corresponding to the desired value of G, and, a(2) is the 
cross-correlation matrix entry the second column of the one by two (1 x 2) matrix 
calculated at Block 208 or Block 209. and stored at Block 210, corresponding to the 
desired value of G. 

[0039] In alternate embodiments, linear interpolation may be utilized to generate 
a particular gamma curve from two or more of the references curves. However, due to 
the non-linear nature of the reference curves, as well as the non-linear nature of the 
particular gamma curve being generated, and the correlation among themselves is also 
non-linear, therefore, linear interpolation is not as accurate as the non-linear 
interpolation techniques described above. 

[0040] Figure 3 is a block diagram of a programmable gamma correction circuit 

300 embodying the principles of the present invention and suitable for such applications 
as gamma correction block 105 of FIGURE 1. 

[0041] In the illustrated embodiment, three reference curves corresponding to 

reference gamma values 1.6, 2.3, and 3.3 are stored in corresponding memory spaces 

301 - 303. The reference gamma values may be stored in memory spaces 301 - 303 
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as LUTs or as piecewise linear approximations of at least 64 segments. Memories 301 
- 303 are indexed by the current value of lumajn. 

[0042] According to the principles of the present invention, gamma curves may 
be generated between a minimum gamma value (e.g. 1 .4) and a maximum gamma 
value (e.g. 3.5). In the preferred embodiment, the locations of the lower and upper 
reference gamma curves are not placed at the exact minima and maxima endpoints, 
and also the middle reference gamma curve is not placed at the exact middle point of 
this range (1.4-3.5) of gamma values. Rather, the lower and upper reference gamma 
curves are placed slightly inside the minimum and maximum endpoints, and the middle 
reference curve location is also slightly different from the exact mid location of the 
above minimum and maximum range. For example, the endpoints of the programmable 
gamma range are at 1 .4 and 3.5 and their midpoint is at 2.45, but the three fixed 
reference curves are placed at gamma values of 1 .6, 3.3 and 2.3, respectively. To 
generate gamma conrection curves with gamma values between 1 .4 and 1 .6, and 3.3 
and 3.5, non-linear extrapolation is performed. Due to the non-linear characteristics of 
the reference curves, the curves being generated, as well as the correlation amongst 
themselves, selecting gamma values for generating the reference curves, which are 
slightly within the maximum and minimum available gamma values and slightly different 
from the mid point of this range, has been found to minimize overall errors. 
[0043] For each value of G programmed in accordance with procedure 200, the 
lower cross-correlation matrix aj is stored in LUT 304, while the corresponding upper 
cross-correlation matrix a_u is stored in LUT 305. 

[0044] A set of selectors (multiplexers) 306 and 307 select, in response to the 
output of comparator 308, upper and lower reference values indexed by the current 
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value of GAMMA G. Specifically, if the current desired input value G, which is stored in 
register 310, is less than or equal to 2.3, then the output from memory 301 is the lower 
reference gamma value and the output from memory 302 is the upper reference gamma 
value. On the other hand, if the current value of G in register 310 Is greater than 2.3, 
then the output from memory 302 becomes the lower reference gamma value and the 
output from memory 303 is the upper reference gamma value. 
[0045] The cunrent value of G in register also indexes the entries of 
corresponding upper and lower cross-correlation matrices aj and a_u stored in LUTs 
304 and 305. The particular matrix accessed is selected by selectors (multiplexers) 31 1 
and 312 under control of the output comparator 313. In particular, when the current 
value of LUMAjN is less than 32, then the indexed matrix from lower LUT 304 is 
accessed, otherwise the indexed matrix from upper LUT 305 is accessed. 
[0046] The gamma con-ected luminance value out is calculated, in accordance 
with Block 216 of procedure 200 of FIGURE 2. by multipliers 314 and 315 and adder 
316. Specifically, for the selected matrix from either LUT 304 or 305, multiplier 314 
multiplies the entry a(1) with the selected lower reference gamma value output from 
multiplexer 306. The access matrix entry a(2) is multiplied by multiplier 315 by the 
upper reference gamma value output from multiplexer 307. The final gamma corrected 
output value out is then generated by summer 316 by taking the sum of the output of 
multipliers 314 and 315. 

[0047] In system 100 shown in FIGURE 1 , the gamma corrected luminance value 
OUT is passed to video interpolation and processing block 106 for further processing, 
[0048] While a particular embodiment of the invention has been shown and 
described, changes and modifications may be made therein without departing from the 
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invention in its broader aspects, and, therefore, the aim in the appended claims is to 
cover all such changes and modifications that fall within the true spirit and scope of the 
invention. 
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